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ABSTRACT 

Is it significant that the intrinsic outputs of several BL Lacs are observed to level off at values of about 10 46 erg s _1 ? In searching for 
an answer, we compare y-ray observations by the AGILE satellite of the BL Lac S5 0716+714 with those of Mrk 421 and Mrk 501; 
the former are particularly marked by intense flares up to fluxes of 2 x 10~ 6 photons cirT 2 s~ ] in the 0.1-10 GeV energy range. These 
"dry" BL Lacs show evidence of neither thermal disk emissions nor emission lines signaling any accreting or surrounding gas; the 
spectral distributions of their pure non-thermal radiations are effectively represented by the synchrotron self-Compton process. With 
source parameters correspondingly derived and tuned with simultaneous multiwavelength observations, we find for S5 0716+714 a 
total jet power of about 3 x 10 45 erg S" 1 , which makes it one of the brightest dry BL Lacs so far detected in y rays. We evaluate the 
mass of the associated Kerr hole to be around 5 x 10 8 Mq, implying that the source is significantly gauged in terms of the maximal 
power around 4 x 10 45 erg s _I extractable via the Blandford-Znajek electrodynamical mechanism; other dry BL Lacs observed in 
y rays remain well below that threshold. These findings and those forthcoming from Fermi-LAT will provide a powerful test of 
electrodynamics in the surroundings of the hole, that are dominated by GR effects. 

Key words. BL Lacertae objects: general - Radiation mechanisms: non-thermal - Black hole physics - Accretion, accretion disks 



1. Introduction 

Blazars rank among the brightest active galactic nuclei on the ba- 
sis of their inferred isotropic luminosities that may attain some 
Lisa ~ 10 48 erg s _1 . Actually, these sources radiate from a nar- 
row relativistic jet closely aligned with the observer's line of 
sight. The jet emits highly beamed non-thermal radiations, with 
observed fluxes enhanced by aberration and Doppler effects of 
' Special Relativity (Begelman, Blandford & Rees I1984I Konigl 
1986; Urry & Padovani 119951 1. So the luminosities L is „ greatly 
exceed the intrinsic outputs of the jets that easily level off at 
10- 2 - 10- 3 L, Jo . 

The BL Lac objects (henceforth BL Lacs) in particular are 
blazars that show no or just weak and intermittent emission lines. 
Their spectra are represented well as a continuous spectral en- 
ergy distribution (SED) S v - vF v featuring two peaks: one at 
a lower frequency due to synchrotron emission by highly rela- 
tivistic electrons; and a higher frequency counterpart due to in- 
verse Compton upscattering by the same electron population of 
seed photons provided by the synchrotron emission itself (syn- 
chrotron self-Compton, SSC; see Jones, O'Dell & Stei n |1974| 
Marscher & Gear [19831 Maraschi, Ghisellini & Celotti [19921 , 
with possible additions from sources external to the beam (ex- 
ternal Compton, EC; see Dermer & Schhckeiser |1993l Sikora, 
Begelman & Rees [T994] l. 

The BL Lacs also exhibit strong variability on timescales of 
days to minutes with substantial flux variations (flares) particu- 
larly at high energies as realized early on (see Setti & Woltjer 
119941 and references therein). 

Here we focus on "dry" BL Lacs, that is, sources with no 
evidence of surrounding gas, such as emission lines or a big 
blue bump (see Peterson 119971 Kembhavi & Narlikar [19991 1 re- 



lated to current accretion. They provide an appropriate testing 
ground for comparing their intrinsic outputs with maximal pow- 
ers extractable from rotating supermassive black holes and from 
the dragged accretion disks by means of large-scale electro- 
magnetic fields, via the intriguing, variously debated Blandford- 
Znajek electrodynamics (BZ, Blandford & Znajek 119771 see 
also Ghosh & Abramowicz 119971 Kr olik|1999l Livio, Ogilvie 
& Pringle [19991 Cava liere & D'Elia 120021 McKinney 120051 
Nemmen et al. 120071 Tchekhovskoy, McKinney & Narayan 
2009). The bare hole contribution can yield up to Lk ~ 2 x 
10 45 (m./io 9 m ) erg s l , given the hole mass M # in units of 

10 9 M Q and a magnetic field B ~ 10 4 G threading its horizon. 

In the following, we adopt the standard, flat cosmology with 
H = 72 km s 1 Mpc" 1 and Q A = 0.74 (Dunkley et al. |25fJ9l . 



2. The radiation process 

The SEDs of the dry BL Lacs are widely understood in terms 
of the simple, homogeneous SSC process. This is based on ra- 
diations produced in a region containing a magnetic field and 
relativistic electrons accelerated to high random energies ymc 2 
(withy up to 10 6 -10 8 )that move toward the observer with bulk 
Lorentz factors r ~ 10 - 20 (see Ghisellini et al. [19931 . 

To begin with, we assume the sources to have an isotropic 
geometry with a radius R as a single size parameter, and to con- 
tain the relativistic electrons and non-relativistic protons with the 
same F, at a common density n. Observed (primed) frequencies 
and fluxes are related to the rest frame (unprimed) quantities by 
means of V — v 6 and F' — F 5 A (Begelman, Blandford & Rees 
1984), where 6 — [T(l -/JcosS)] 1 is the beaming factor re- 
lated to the angle 6 between the jet and the line of sight. Small 
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viewing angles G ~ l/T yield 6 « 2F. Correspondingly, the in- 
trinsic variability and crossing times R/c will be longer than the 
observed ones R/6c. 

On empirical and theoretical grounds, we adopt log- 
parabolic shapes for the electron energy distributions. These are 
obtained from a Fokker-Planck equation in the presence of sys- 
tematic and stochastic acceleration processes as first shown by 
Kardashev (1 19621 1 and computed in detail by Paggi et al. (2009); 
the acceleration times scale as t a <x y/E, in terms of the effective 
electric field E (Cavaliere & D'Elia 2002). We therefore write 
the electron distribution in the form 
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N U 



-s-r log 



(1) 



where s is the constant contribution to the slope, r is the "cur- 
vature", and 7q is the injection energy. The energetic content of 
such an electron population can be expressed in terms of the rms 

adimensional energy y p = y 2 N(y) dy/ fN(y)dy) , to yield 

an energy density close to n y p mc 2 . 

The emitted synchrotron SED is correspondingly given by 
(Massaro et al. 120041 Tramacere et al. 120071 and references 
therein) 



Sy - Sf 
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with a constant contribution a - (s - 3)/2 to the spectral index, a 
spectral curvature b « r/5, and a peak frequency v$ oc By 2 . For 
the IC, an analogous SED applies, where in the Thomson regime 
one has a = (s - 3)/2, b « r/10 and a peak frequency vc oc By A p , 
whilst in the Klein-Nishina regime a — s, b « r, and vq oc y p 
hold. 

In these SEDs and distributions, the peaks where most of 
the energy resides are not materially affected by radiative cool- 
ing active on timescales t c oc l/y p longer than the crossing time 
R/c; cooling will rather erode the high energy tails. As an added 
bonus, the synchrotron and IC radiations from a log-parabolic 
electron population irreversibly broaden under the action of the 
stochastic acceleration component, following b oc l/t (Paggi et 
al.|2009); thus, a sudden increase in the spectral curvature will 
mark the emergence of a new electron population. Last but not 
least, Eq. |2]closely fits (as we illustrate in Fig. |T} the spectra of 
the sources that we focus on in view of extended spectral cov- 
erage provided by their multiwavelength, simultaneous observa- 
tions. 

This is the case for S5 0716+714, which has the follow- 
ing data available: a low state provided by historical AIT and 
EGRET data (Lin et al. [19951 ); a high state in September 2007, 
in the y-ray range covered by AGILE-GRID where the IC peaks 
for this "intermediate" BL Lac (Padovani & Giommi 1995); op- 
tical and radio data taken with GASP-WEBT on September 7-12 
(Villata et al. 120081 Vittorini et al. 120091). Th e multiwavelength 
variations observed by Giommi et al. (2008 ) and the increased 
spectral curvature (see Fig.[T) in 2007 are indicative of the injec- 
tion of a second component. 

The source may be compared with Mrk 501, for which we 
consider two states with simultaneous BeppoSAX and CAT ob- 
servations on April 7 and 16, 1997 (Massaro et al. 2006); and 
with Mrk 421, for which we have low and high states in 2000 
from BeppoSAX and HEGRA data (Konopelko et al. 120031 ), and 
multiwavelength observations performed with GASP-WEBT, 
RXTE/ASM, Swift, SuperAGILE, AGILE-GRID, ARGO-YBJ, 
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Fig. 1. SEDs of the BL Lacs considered in the text: S5 0716+174 
(upper-left frame), Mrk 501 (upper-right frame) and Mrk 421 
(lower-left and lower-right frames), each in a low (blue line) 
and a high (red line) state. In terms of intrinsic luminosity S5 
0716+714 is the strongest source (data referenced in Sect. [2]). 



and VERITAS in June 2008 (Donnarumma et al. [20091 Di 
Sciascio et al. 120091 ) . 



3. The source power 

We are interested in the intrinsic outputs referred to the jet frame, 
rather than in the luminosities L, JO = AnD 2 L F inferred from in- 
sisting on an isotropic distribution of the observed flux F, at the 
luminosity distance Di. We assume one "cold" proton per elec- 
tron satisfying (y) £ m p /m e (with the average (y) — y p x 10~^ 4r 
bounded in terms of the electron m e and the proton m p masses), 
and follow Celotti & Ghisellini (2008 ) in writing for the intrinsic 
radiative luminosity L r contributed by both the synchrotron and 
IC radiations and for the related powers carried by the jet, the 
expressions 



L r = L iso T 2 /6 4 ~ L is0 /16T 2 , 
L e = —nR 2 cnm e c 2 (y)Y 2 , 
L p ~ L e m p /m e (y), L B < 



(3) 
(4) 
(5) 



The total jet power is therefore given by Lj - L r + L e + L p + Lb\ 
this is dominated by L r and by L e , with Lb < L p % L e for the 
fields B ~ 0.1 — 1 G implied by the spectral fits. 

The simultaneous, multiwavelength observations enable ex- 
tended spectral fits to determine the five key observables (beside 
the spectral curvature b) from the SSC model, namely: the syn- 
chrotron peak frequency and flux, the IC peak frequency and 
flux, and the variation time (see Paggi et al. 2009 ). These lead to 
robust evaluations of the five source parameters n, R, B, (y), and 
6 (or T) entering Eqs.[3]-|5j the main parameters are collected in 
TableQ] 



4. The BZ benchmark 

As anticipated in Sect. [1] a natural benchmark for these pow- 
ers is provided by the BZ mechanism for electrodynamical en- 
ergy extraction from a Kerr hole spun up to maximal rotation by 
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Source name 


S 


y p 


L r 


L T 




S5 0716+714 (low) 


12 


9.1 x 10 2 


4.2 x 10 44 


1.3 x 10 45 




S5 0716+714 (high) 


19 


2.8 x 10 3 


1.5 x 10 45 


3.1 x 10 45 


10-, 


Mrk 501 (low) 


15 


1.4 x 10 5 


4.9 x 10 42 


1.2 x 10 43 




Mrk 501 (high) 


15 


1.9 x 10 5 


2.3 x 10 43 


3.1 x 10 43 


1 -.- 


Mrk 421 (2000, low) 


20 


5.2 x 10 4 


4.0 x 10 42 


1.9 x 10 43 




Mrk 421 (2000, high) 


20 


6.1 x 10 4 


8.1 x 10 42 


3.1 x 10 43 




Mrk 421 (2008, low) 


20 


2.5 x 10 4 


1.6 x 10 43 


4.1 x 10 43 


0.1 - 


Mrk 421 (2008, high) 


20 


3.4 x 10 4 


2.1 x 10 43 


4.6 x 10 43 


si : 
m : 



Table 1. Parameters for the 
The redshifts are z — 0.03 1 
501 and z = 0.030 for Mrk 



BL Lac sources discussed in the text, 
for S5 0716+714, z = 0.034 for Mrk 
412; L r and Lj are given in erg s _1 . 



past accretion episodes. A minimal, vestigial disk is required to 
hold the poloidal magnetic field threading the horizon; the disk is 
kept active by low accretion rates m ^ 10~ 2 in Eddington units, 
loses angular momentum mainly via the large-scale field, and 
contributes some 3 L% to the total power (Blandford & Znajek 
119771 Livio et al. H999t . The two contributions add to yield 



L BZ * 8 x 10 4 



10 9 M P 



ergs 



(6) 



We note that the balance B 2 /4n ~ p between the magnetic stress 
and the kinetic or radiation pressure p in the disk yields B ~ 
10 4 G; for a radiation-pressure dominated disk, we have at the 
inner rim B 2 oc 1 /M # , so B has dropped out of Eq. [6] 

The hole mass is then the key parameter, that we evaluate 
from its correlation with the absolute red magnitude Mr of the 
host galactic bulge (Ferrarese et al. 120001 Gebhardt et al. 2000; 
Falomo et al. 120031 : for our cosmology this reads 



log (^r- a50M *- 2 - 61 ' 



(7) 



with scatter +0.4 dex (Bettoni et al. 12003 ). For the host galaxy 
of S5 0716+714, observat ions o f the magnitude R = 18.3 + 0.5 
reported by Nilsson et al. d2008l l, besides indicating the redshift 
z = 0.31 + 0.08, yield a mass M. 5.5+f-° 10 s M Q ; the central 
value is consistent with estimates from microvariability of the 
optical flux (Sasada et al. 120081 , which yield M. ~ 10 8 M o . For 
Mrk 501 and Mrk 421, one obtains M. 1.0+J 4 , 10 9 M o and 
M. a 10 s M Q , respectively. 

Our results normalized to the respective Lbz from Eq.|6]are 
represented in Fig. [2] During flares, the electron rms energies 
(and the peak frequencies) are boosted in all sources, and so are 
the luminosities; this indicates that rising flares are directly re- 
lated to increased acceleration of the emitting electrons. 

We emphasize that the powerful source S5 0716+714 is ap- 
parently constrained to move sideways, as if to skirt the BZ limit; 
Mrk 421 in the 2008 states exhibits a similar behaviour, although 
with lower significance. On the other hand, the weaker source 
Mrk 421 in 2000 and the yet weaker Mrk 501 remain consider- 
ably below the BZ limit and so are expected to be free to move 
more vertically, as they do. 



5. Discussion and conclusions 



For dry BL Lacs with accretion rates m < 10~ 2 , the SSC radia- 
tion process provides a robust evaluation of the jet luminosities. 
Whence we conclude that Lbz provides a significant benchmark 
for the output of the BL Lacs discussed here, and indeed an up- 
per limit to both their quiescent states and flares. In fact, during 




Mrk 421 (08) 




Mrk 421 (00) 



1 



Fig. 2. For the named sources the total jet luminosities normal- 
ized to their BZ power are plotted against the electron rms en- 
ergy. Bars represent the hole mass uncertainties reflecting those 
(at the l—o" level) in host bulge luminosities, and including the 
scatter in Eq. (0. 



a recent flare S5 0716+714 was observed to be constrained by 
Lbz 55 10 46 erg s _1 , and a similar behavior was observed in 2008 
for Mrk 42 1 . Non-thermal, beamed powers in the range Lk -Lbz 
also provide evidence of an accretion disk that is active mainly in 
launching and channeling the jets by means of large-scale fields. 

Referring to Fig.[3]and its caption, we note that during flares 
the sources move in the Lj - y pea k plane away from the envelope 
that is outlined by bright BL Lacs with increasing rates m; the 
envelope ends up in the locus of the yet brighter Flat Spectrum 
Radio Quasars (FSRQs) with m ~ 1 . The flares then move into a 
region of faster radiative cooling (Celotti & Ghisellini 2008J and 
references therein). This implies short-lived flares on timescales 
S5 1 day, or requires shorter acceleration times t a ~ y/E with 
higher E, as an alternative to structured sources such as de- 
celerating (Georganopoulos & Kazanas 120031) or spine-sheath 
jets with inner scale Ri < R (Tavecchio & Ghiselllini 2008). 
Faster acceleration and deviations from the envelope are consis- 
tent with flares caused by electron boost rather than episodes of 
increased accretion onto the disk. 

In this context, we recall (see Blandford 1990; Padovani et 
al. 2007 ) that sources lying along the envelope in Fig.[3]at higher 
L and lower y pea k often exhibit stronger evidence of current ac- 
cretion up to m ~ 1, such as thermal emissions and surrounding 
gas (big blue bump and broad emission lines), with a larger con- 
tribution from EC. In fact, the progression from dry BL Lacs 
to FSRQs is likely to involve an enhanced and extended disk 
contribution as described by Blandford & Payne (1982), start- 
ing with "wet" BL Lacs with m ~ 10~'; these feature larger EC 
contributions (Dermer et al. 2008 ) and looming evidence of gas, 
including some thermal disk emission and weak or intermittent 
lines (Celotti, Ghisellini & Fabian 2007). The last step in this 
progression is constituted by the powerful FSRQs with extant 
broad lines, a big blue bump from disks accreting at full rates 
m ~ 1 , and a dominant or towering EC (Maraschi & Tavecchio 
I200TT) . 

We add that the outputs of even misaligned BL Lacs may 
be calorimetrically gauged from their feedback actions on the 
intra-cluster plasma surrounding their host galaxy when located 
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Fig. 3. Bright BL Lacs in their context, adapted from Celotti & Ghisellini (120081) with historical data in terms of Lj and the energy 
J peak related to the synchrotron peak (y pea k = J P X lO 1 ^ 2 '). Blue circles indicate dry, while violet circles indicate wet BL Lacs. The 
lower-left region of the diagram corresponds to the source condition t c > R/c, and the upper-right to t c > t a . Bright FSRQs lie at 
lower jpeak and higher powers, with increasing signs of current m — » 1; selection effects depopulate weaker sources in the lower 
region (see Padovani |2007l l. The dashed oval highlights sources in the transition region from dry to wet BL Lacs, interesting to 
compare with Lgz from Eq.|6]particularly when at z > 0.3. 



in a cluster or a group of galaxies, as discussed by McNamara 
et al. ( 120091 1. These authors evaluate average powers around 
10 46 erg s" 1 injected into the cluster MS0735. 6+7421, and pos- 
sibly also in the cluster A2029 and the group AWM 4. 

The whole of the above evidence provides observational sup- 
port to the relevance of the electrodynamical BZ mechanism, 
and invites extended sampling of other interesting sources (see 
Fig.E). 

If in dry BL Lacs with M, < 10 9 M© the L B z limit were 
found to be substantially exceeded by outputs Lj > 10 45 erg s _1 , 
this would require B > 10 4 G at the Kerr horizon. These fields 
imply large dynamical stresses bounded only by B 2 jAn < pc 2 , 
associated with particle orbits plunging from the disk toward 
the hole horizon (Meier 2002]) into a region fully controlled by 
strong gravity effects. 

Thus, all such sources will provide powerful tests for the 
coupling of electrodynamics with General Relativity in full 
swing, and constitute an exciting arena for AGILE and Fermi- 
LAT data. 
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